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Q^' Abstract 

«!-> : 

,^ , The gravitational dual of an insulator/superconductor transition driven by increas- 

ing the chemical potential has recently been constructed. However, the system was 

CN ' studied in a probe limit and only a part of the phase diagram was obtained. We in- 

clude the backreaction and construct the complete phase diagram for this system. For 
fixed chemical potential there are typically two phase transitions as the temperature 

r^ ' is lowered. Surprisingly, for a certain range of parameters, the system first becomes a 

superconductor and then becomes an insulator as the temperature approaches zero. As 

LT* I a byproduct of our analysis, we also construct the gravitational dual of a Bose-Einstein 

f— ^ ' condensate of glueballs in a confining gauge theory. 
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1 Introduction 

The AdS/CFT correspondence [HE], [3] has provided a framework where strongly coupled field 
theories can be studied using a weakly coupled gravity dual. Recently, this correspondence 
has led to gravitational descriptions of superfiuids and superconductors [H El E]- In the 
simplest realization, a black hole in an Einstein-Maxwell-charged scalar theory becomes 
unstable to forming scalar hair below a critical temperature T^. In the dual field theory, this 
instability corresponds to a second order phase transition which spontaneously breaks a U{1) 
symmetry. In this way, the field theory exhibits many properties shared by superconductors 
and superfiuids. This basic construction has been extended in various ways (for reviews see [71 
[8l[9l[T0]). Such systems are called holographic superconductors. Above the phase transition, 
the black hole describes a conductor, so this is a standard conductor/superconductor phase 
transition. 

In a recent paper [11], it was shown that the same simple bulk theory with modified 
boundary conditions can describe an insulator/superconductor transition. The new bound- 
ary conditions are simply that one direction is compactified to a circle (and fermions are 
antiperiodic around this circle). The insulator is described by the AdS soliton [T3l [H], 
which has a mass gap and has previously been used to describe a confining vacuum of a dual 
gauge theory. In the well studied example of asymptotically AdS^ x S^ solutions which are 
dual to A/" = 4 super Yang-Mills, the new boundary conditions give all fermions masses at 
tree level and the scalars pick up masses at one loop. The low energy theory is therefore a 
pure 2 + 1-dimensional gauge theory which has a confining vacuum. The AdS soliton is the 
gravitational description of this state. For our purposes, we will not consider any particular 
embedding in string theory and instead just view the soliton as describing an insulator. 

It was shown in [TT] that if one adds a chemical potential /i to the AdS soliton, there 
is a second order phase transition at a critical value fie beyond which the charged scalar 
field turns on, even at zero temperature. It was also shown that for all fi > fie, the DC 
conductivity becomes infinite. So this is an insulator/superconductor transition. 

For small fi, it is known that there is another phase transition when one increases the 
temperature [12] ■ This is a first order phase transition between the AdS soliton and the 
AdS black hole). In gauge theory language, this is a confining/deconfining phase transition 
|13] . while in the condensed matter interpretation it is a insulator/conductor transition. 
Therefore, there are a total of four phases in this system given by the AdS soliton, the AdS 



^This phase transition is the planar analogue of the Hawking-Page transition for global AdS [15] 



black hole, and their superconducting phases. Thus, just by compactifying one direction, 
the simple holographic superconductor model gets a much richer phase diagram. 

The analysis in [H] was done in the probe limit where the backreaction of the matter 
fields is ignored. This is justified in the limit of large charge q on the scalar field with g/x 
fixed. The transition in which the soliton develops a nonzero scalar field was studied in this 
limit, but since the probe limit requires /i ^ 1, the full phase diagram was not constructed. 

In this paper, we extend this previous work and include the backreaction of the matter 
fields on the soliton geometry and black hole. We also compute the full phase diagrams as 
a function of temperature and chemical potential for various q. These are given in Figure 6. 
For all g > 1, we find that the black hole phases never dominate at low temperature. The 
system always prefers either the soliton insulator or soliton superconductor as T — )■ 0. The 
four phases typically meet in two triple points where three of the phases can coexist. As we 
lower g, the two triple points merge into one "quadruple" point and then separate, passing 
through each other. This produces a new phase boundary in which the system goes from a 
conductor to a superconductor and then to an insulator as the temperature is lowered! This 
is very different from the superconductor/insulator transition studied in |Tl]. In the previous 
case, the transition was second order and occurred by decreasing the chemical potential at 
fixed temperature. The new superconductor/insulator transition is first order and occurs 
when decreasing the temperature at fixed chemical potential. 

As an aside, we also consider the limiting case of a neutral scalar field. When g = 0, the 
AdS soliton is never unstable to turning on a static scalar field. However, there are finite 
energy solutions in which the geometry is static and the scalar field has pure harmonic time 
dependence. These are planar analogues of boson stars. (For a review of boson stars see [16j . 
and for a discussion of spherical boson stars in asymptotically AdS spacetimes see [I7j.) Like 
typical boson star solutions, we find that there is a maximum mass that can be supported. 
In terms of a dual gauge theory, these solutions correspond to Bose-Einstein condensates of 
glueballs. 

The paper is organized as follows. In the following section, we write down the equations 
of motion for our systems and discuss the solutions with no scalar field. In section three, 
we numerically construct the soliton solutions with scalar backreaction, including the g = 
case. We then construct the black hole solutions with scalar hair in section four. In section 
five, we assemble the phase diagrams of our system. The last section summarizes our results. 



2 The Bulk Equations and Solutions without Scalar 

We begin with a five-dimensional Einstein-Maxwell-scalar action: 



s = I d'xy^ ( ^ + 7^ - l^'^'^F^'^ - w\' - ^vr 



(2.1) 



Here, we are writing F = dA and D^ = V^ — iqA^. The equations of motion obtained from 
this action are the scalar equation 
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Maxwell's equations 
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and Einstein's equations 
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First, let's consider solutions with ip = Q and At = 0(r). One solution with planar symmetry 
is the AdS Reissner- Nordstrom black hole, given by [7J 
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We have written the metric in this form so that r+ is the horizon radius and fi is the chemical 
potential of the dual field theory. The temperature of this black hole is 
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We will later want to make z a periodic coordinate. Any period can be chosen for z. 



There is another planar solution, the AdS soliton, which is given by 

ds^ = ^L^ + r^dx" + dy^ - de) + f{r)dii\ f{r) = r^ - 4 ■ (2-9) 



fir) ' r 

and = 0. This solution can be obtained by a double Wick rotation of the AdS Schwarzschild 
black hole. That is, we set /i = in (12. 5p and (12. 6 p and make the substitutions t ^ irj and 
z — !• it. The geometry resembles a cigar with the tip at r = tq. Unlike the black hole, in 
order to avoid a conical singularity at the tip, r] must be chosen to have a period of 

7 = - . (2.10) 

To 

The soliton can be considered at any temperature. That is, we are free to give imaginary 
time any period /3 = 1/T we wish. 

Due to the periodicity of r], the theory dual to the AdS soliton lives on a space with 
one direction compactified to a circle with length 7. Because the AdS soliton only exists for 
r > To, this dual field theory is confining and has a mass gap. Consider this theory at finite 
temperature. There are two phases given by the Euclidean black hole (with fi = and z 
periodic with period 7) and the Euclidean soliton. These are actually the same geometry, 
differing only in the interpretation of which S^ at infinity is Euclidean time and which is a 
spacelike direction. For the black hole, the Euclidean time direction is contractible, while 
for the soliton, the spatial direction is contractible. Since the free energy is computed from 
the Euclidean action, there is a phase transition when the two circles have the same length: 
/3 = 7. This is a first order transition, exactly analogous to the Hawking-Page transition for 
global AdS [T5]. 

One can extend this to nonzero chemical potential by noticing that one can add (f) = fi 
to the soliton solution without changing the metric. As we increase fi, there is an instability 
which turns on the charged scalar field (in both the soliton and black hole phases). This 
was analyzed in the probe limit (g — )■ 00, keeping qip and q(f) fixed) in which backreaction on 
the metric is ignored in flT]. We would like to move beyond the probe limit and study the 
phase transitions with full backreaction. To do this, we must solve the full set of differential 
equations, which we turn to in the following sections. 



3 AdS Soliton with Scalar 

3.1 Equations of Motion and Boundary Conditions 

Now let us now consider the soliton solutions with ^ 7^ 0. We will work in units in which 
L = 1. We are interested in including backreaction so we choose the metric ansatz 

ds^ = r^ (e^(^-^B{r)dr]^ + dx^ + dy^ - e^^cit^) + -^ . (3.1) 

\ / r^i5(r) 

We require that B vanish at some radius tq, which is the tip of the soliton. As in the case 
without a scalar field, smoothness at the tip requires that r] be periodic with period 

We would like to consider solutions of the form 

At = 0(r), ^P = ilj{r)e-"^' . (3.3) 

For q ^ 0, there is a symmetry 

^^V'e'"*, 0^0 + -, (3.4) 

Q 

that lets us set u = 0. If g = 0, then decouples from ip and we cannot gauge away u. 
With this ansatz, the scalar and Maxwell equations become 

^'+- + — + ^ + — U'+^TT —^r^-^^ - m^ U = , (3.5) 
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For simplicity, we choose m? = — ^. Our analysis can be generalized to any mass above 
the Breitenlohner-Freedman bound [18] m^ > — ~^ = —4. Near the boundary r = 00, 
the scalar and Maxwell equations have the form 

ijrf-J/ ii 'Y' / 



= /i-^ + .... (3.8) 

The quantities /x and p are interpreted as the chemical potential and the charge density in 
the dual field theory, respectively. At our value of m?^ the V'*^^-' and ■i/'^^^ terms are both 
normalizable so we have a choice of boundary conditions. For numerical simplicity, we only 
consider the case i})'^^'^ = 0. The constants ip^^^ and ^^"^^ can be used to define operators on 
the dual field theory with mass dimension A = 3/2 and A = 5/2, respectively. Up to a 
normalization, the expectation value of the scalar operator O is given by 

(C) = ^(2) . (3.9) 

The nontrivial components of Einstein's equations are the tt, rr, rji], and xx. These 
equations are not all independent. The xx component, for example, can be derived from 
combinations of the other equation of motion and their derivatives. Therefore, we only need 
three linearly independent combinations of these equations. Subtracting the rji] equation 
from the rr equation allows us to solve for A': 

A' = ^, —- . (3.10) 

r(6 + rC) ^ ' 

The remaining two equations come from the xx — tt component 

''" - ¥" - (; - T - f ) ^' - (^" - ^^^^1^) ^ - « ("^) 

and rjrj + tt — xx 

(3.12) 
Note that we can substitute (I3.10p into the other equations of motion to eliminate A. This 
does not change the order of our differential equations. We can then integrate (I3.10p sepa- 
rately after solving the other equations. 

By considering a series solution about the tip of the soliton, r = ro, and using -B(ro) = 0, 
we are left with four independent parameters: ro, ^(^o), 0(?"o)5 ^'('"o)- However, there are 
two useful symmetries. The scaling symmetry 

r ^ ar, {ri,x,y,t) ^ {ri,x,y,t)/a, (j) ^ acj) (3.13) 



can be used to set tq = 1 for performing numerics. In general, different solutions obtained 
this way will have different periods 7 for the r] coordinate. In order to make appropriate com- 
parisons between different solutions, the boundary geometry must be the same. Therefore, 
we will later use f l3.13p again to set all of the periods 7 equal. Similarly, the symmetry 

e"^ -^ a^e"^, t ^ at, (/> -> (j)/a (3.14) 

lets us pick an arbitrary value of C(ro). After solving the differential equations, we can read 
off the value of C(oo) and use this symmetry to satisfy the asymptotic condition C(oo) = 0. 

If g 7^ 0, we can choose 0(ro) as a shooting parameter, which we tune to satisfy the 
remaining asymptotic condition, ip^^^ = 0. li q = 0, (j) decouples from ip (and will be set to 
zero). However precisely in this case, u cannot be gauged away and becomes our shooting 
parameter. 

So for each choice of q and ipi'^^o), we can solve the equations of motion for ip, 0, B, and C. 
The remaining function A can be found by integrating (I3.10p with the condition y4(oo) = 0. 
Thus for fixed q, we have a one-parameter family of solutions. After solving the equations, 
the AdS/CFT correspondence allows us to obtain the chemical potential /i and the charge 
density p from the asymptotic behaviour of (p according to (13.81) . We can also obtain (O) 
through the asymptotic behaviour of ip through (13. 7p and (13. 9p . 

3.2 q^O 

Our main interest is the case when the scalar field is charged: q ^ 0- One finds that when fi 
is small enough, there are no solutions for ip which satisfy our boundary condition. However 
when one increases /i, such solutions exist. It is convenient to express the results in terms 
of scale-invariant quantities. The condensate (O) has dimension 5/2, the chemical potential 
has dimension one, and the charge density has dimension three. Using the length of the 
circle at infinity, 7, as a scale we plot ((9)^/^7 and pj^ as functions of p'j in Figure [1] for 
q = 2. We have found that the plots for q = 10, 5, and 1.2 are qualitatively similar. 

The main conclusion of Figure 1 is that a phase transition occurs at a critical value of 
p. The transition is second-order, which can be seen from the p-p plot. If one computes the 
optical conductivity in the AdS soliton background (with no scalar) one finds that Re[(j(a;)] 
is a series of delta functions at nonzero frequency [H]. In particular, the DC conductivity 
is zero, so this describes an insulator. When the scalar field is nonzero, there is an ad- 
ditional delta function at a; = indicating the onset of superconductivity pjj. Thus, we 
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Figure 1: The value of the condensate (left), and the charge density (right) as a function of chemical 
potential, for the soliton with q = 2. 



have a holographic superconductor/insulator phase transition that can exist at strictly zero 
temperature. 

There is an additional complication when q < 1.2. Recall that our one parameter family of 
solutions can be labelled by the value of the scalar field at the tip of the soliton. For q > 1.2, 
as one increases '?/'(ro), /i7 monotonically increases. So every value of the chemical potential 
has at most one soliton solution with scalar. When q < 1.2, fi'j does not monotonically 
increase with ip{ro) and there are up to three different solutions for the same chemical 
potential. The free energjo versus chemical potential plots develop "swallow tails" (see 
Figure [2]) common in first order phase transitions. A similar effect has been observed in 
spherical Reissner- Nordstrom AdS black holes (when the total charge is fixed rather than the 
chemical potential) [I9]. This new phase transition first appears inside the superconducting 
phase at some fidis- The charge density p and condensate (O) are also discontinuous at 
fidis (see Figure [3]). It is not clear what physically distinguishes these two phases of the 
superconductor. As one lowers q, the first order phase transition moves over and coincides 
with the superconductor/insulator transition when q ^ 1. 



3.3 



q 







Now let's consider g = and turn off the chemical potential: fi = 0. The Maxwell field 
drops out of the problem and we are left with just gravity coupled to a complex scalar field. 
We look for solutions of the form ip = il){r)e~'^^^ with cj 7^ 0. For every choice of V'(^o); 
one can adjust uj so that ?/'(r) is normalizable at infinity. Therefore, our solutions are boson 
stars in the sense that we have a static spacetime solution with a complex scalar varying 



^See section 5 for more discussion on the computation of free energy. 






Figure 2: The free energy densities of the sohton with scalar (sohd purple), and the soliton without 
scalar (dashed blue), against chemical potential with q = 1.2 (top left), q = 1.15 (top 
right), q = 1.1 (bottom left), and q = 1 (bottom right). Here, we have scaled 7 = tt. 
The q = 1.2 plot shows the typical secord-order phase transition seen in the probe limit. 
The q = 1.15 and q = 1.1 plots have second order phase transitions, but a discontinuity 
within the superconducting phase. The q = 1 plot shows a first order phase transition. 

harmonically with time. Our solutions, however, have planar symmetry while most boson 
stars that have been studied have spherical symmetry. 

The energy of the boson stars are plotted in Figure HI The AdS soliton itself has negative 
energy, i.e., energy below AdS^. As we increase ip{rQ) the energy grows to a positive maximum 
and then decreases. This is similar to previous studies of boson star solutions, in that there 
is a maximum mass that can be supported. We find that the condensate also has a maximum 
value. These results are shown in Figure HI 

In typical boson stars, ipc, the value of ipi'^^o) corresponding to the maximum mass, is 
a boundary between stable and unstable solutions. Solutions with ip{ro) < ipc are stable 
while solutions with ^(rg) > tpc sue unstable [TTj. The unstable solutions will evolve into a 
Schwarzschild black hole. We believe that our solutions will behave similarly, although we 
do not analyze stability here. One simple check is that the energy in the unstable branch is 
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Figure 3: The value of the condensate (left plots), and the charge density (right plots) as a function 
of chemical potential, for the AdS soliton with q = 1.15 (top plots) and g = 1 (bottom 
plots). Here, we have scaled 7 = vr. Note that the q = 1.15 plots have a second order phase 
transition into the superconducting phase (at ^ ^ 1-71), but within the superconducting 
phase (at /i ~ 1-8), there's a discontinuity in the charged density and condensate. The 
q = 1 plots give a first order phase transition. 

always positive and hence there is a black hole with the same energy for it to evolve into. 

When ?/'(ro) is small, the backreaction is negligible and one is solving the linearized 
scalar equations on the AdS soliton background. In terms of the dual gauge theory, this 
describes a spin glueball state above the confining vacuum. By increasing ?/'(ro) and 
including backreaction, one is describing the gravitational dual of a Bose- Einstein condensate 
of glueballs. Figure 4 shows that these condensates remain coherent and do not thermalize 
even when their energy is nearly double that required to drive a deconfinement transition. 

Since the double analytic continuation of the AdS soliton is a planar black hole, it is 
natural to ask what happens if one does the same analytic continuation on the boson star. 
Setting t -^ iz and rj — )■ it, the scalar becomes ilj{r, z) = il){r)e^^ and the metric becomes 
a translationally invariant black hole. It thus appears that one can have inhomogeneous 
scalar hair outside a homogeneous black hole! To understand this strange solution, we must 
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take a closer look at the behavior of the Lagrangian under analytic continuation. The scalar 
field in the boson star is complex, so the analytic continuation results in two real scalar 
fields. Starting with ipi = %p{r) cosut and '?/^2 = V'(?") sinut we see that ipi is the expected 
exponentially growing solution outside the horizon, but tp2 = iipir) sinhuz is pure imaginary. 
This means that it has negative kinetic energy and the theory does not have a stable ground 
state. 




l*(ro) 




lAdo) 



Figure 4: The energy density (left), and the condensate (right) as a function of scalar field at the 
tip, for the AdS soliton with neutral scalar field excited. Here, we have scaled the period 
of the r] coordinate to be 7 = vr. These solutions have a maximum mass and condensate. 



4 AdS Charged Black Hole with Scalar 

In this section we construct asymptotically AdS charged black holes with scalar hair in five 
dimensions with m'^ = —15/4 and various values of q. It is convenient to use the metric 
ansatz 

(4.1) 



ds'' = -g{r)e-''^'Mt' + — - + r\dx^ + dy^ + dz^) 

9{r) 



which has a temperature 



and entropy 



T 



^'(r+)e-x('-+)/2 
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S = AnrlVs = A'nr\-iV2 



(4.2) 



(4.3) 



where the horizon, r+, is defined by (7(r+) = 0. To match with the soliton, we will periodically 
identify the z coordinate with period 7 = vr. Then with the same ansatz for A^ and ip as 
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the soliton (13.3p . and setting u = 0, the scalar and Maxwell's equations are 

V/'+fi-^ + £:),.' + i('f!«!^-,„A,. = o. (4.4) 

\r 2 g J g \ g J 

while the tt and tt — rr components of Einstein's equations become 



^^ 3? 



X- + ^ + ^'\%^^ =0, (4.6) 



The xx component of Einstein's equation can be derived by differentiating the other equa- 
tions of motion. As before, we choose rn? = — ^- Note that (13. 7p and (13.81) still hold in this 
case. 

Now we consider boundary conditions. In addition to g{r+) = 0, we must impose 0(r+) = 
in order for g^^ A^A^ to remain finite at the horizon. Then the independent parameters at 
the horizon are r+, '?/'(r+), 0'(r+), and x(^+)- As in the case for the soliton, we can use the 
scaling symmetries 

r — )■ ar, (r, x, y, t) — )■ (r, x, y, t)/a, g — )■ a^g, (p ^ acj) , (4.8) 

and 

e^ -^ a^e^, t -^ at, (j) -^ (f)/a (4.9) 

to set r+ = 1 and x(oo) = 0. We will later use ( 14. 8 p again to rescale fi in order to make 
comparisons with the soliton solutions. We choose 0'(r+) as a shooting parameter to set 
t/^W = 0. 

After solving the differential equations numerically, we can plot the scale-invariant quan- 
tities (C)^/^//i and p/n^ versus T/jj, as shown in Figure |5l These scale invariant quantities 
are equivalent to those scaled to /i = 1. These results are similar to results for the four- 
dimensional superconductors studied in [5l E]. The scalar turns on quickly below some 
critical temperature through a second-order phase transition. Unlike the AdS soliton, the 
phase transition remains second order even for g < 1. (We have checked as low as q = 0.5). 



12 



<0>^~l'l^ 4 - 































/: 












-..._^ 




1.35 








~— 














"\ 


1.30 














/ : 


- 












\ 


1.25 












/ 


: 


- 














, 3 1.20 
1.15 





__ 








/ 


: 


- 














1.10 
1.05 






^ 


X 


/ 




: 


- 














1.00 








\ 


y 






0.00 


0.02 


0.04 


0.06 0.08 


0.10 


0.12 


0.14 




0.08 


0.10 


0.12 


0.14 
T/;i 


0.16 


0.18 


0.2C 



Figure 5: The value of the condensate (left), and the charge density (right) as a function of chemical 
potential, for the AdS black hole with q = 2. 

5 Phase Diagrams 

We have seen that the system dual to our bulk Einstein-Maxwell-scalar theory has four 
phases. There is an insulating phase described by the AdS soliton, a conducting phase de- 
scribed by the charged black hole, a superconducting phase described by the hairy black 
hole, and a second superconducting phase described by the soliton with nonzero scalar field. 
With slight abuse of terminology, we will refer to the last two phases as the black hole su- 
perconductor and soliton superconductor respectively. A key difference between these two 
superconducting phases is that for the black hole superconductor, the low temperature op- 
tical conductivity has a gap at low frequency but approaches the normal state (nonzero) 
conductivity at larger frequency. For the soliton superconductor, the real part of the con- 
ductivity is just a series of delta functions. 

In this section, we construct the complete phase diagram as a function of temperature 
and chemical potential. There is a third scale in the problem, the size of the compact 
direction at infinity, which we fix to be 7 = vr. This removes the usual scaling symmetry. 
To construct the phase diagram, we must compute the Gibbs free energy of each system in 
the grand canonical ensemble. This can be obtained from the classical gravity action (12.11) 
of the Euclidean solutions: 



lG = m = f3{E -TS- fiQ) 



(5.1) 



where E is the total energy. Written in terms of densities. 
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(5.2) 
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Let's first consider the AdS soliton without the scalar field. The entropy and charge 
density vanish for this solution so the free energy is given (in our units) bjo 

Fsi = esi = -^ = -1 , (5.3) 

7 

where in the last equality we have used the fact that 7 = tt. The free energy of the AdS 
soliton with a charge scalar is given by 



slsc 



esisc - fJ'Psisc ■ (5.4) 



F>^h = -rl{l + ^] ■ (5.5) 



The free energy of these solutions is independent of temperature, so a phase boundary 
between them appears along a line of constant chemical potential. 

Now consider the Reissner-Nordstrom AdS black hole. We can calculate the free energy 
through 

3r2 

We can solve (12. 8p for r+ and substitute the result into (15. 5p to obtain F^h as a function of 
temperature and chemical potential. With the scalar hair turned on, the free energy of the 
black hole is given by 

Fbhsc = ^bhsc — Ts — fiPbhsc , (5.6) 

where the temperature and entropy can be obtained from (14.21) and (14. 3p . Since the critical 
temperature for the black hole to develop scalar hair is proportional to /i, this phase boundary 
is a line that passes through the point {p,T) = (0,0). 

Now we can begin to assemble our T-fi phase diagram. The lowest value of F among our 
four phases is favored. We assemble the phase diagram by the following procedure. First, 
we analyze our system at a fixed value of /i. We must use (14.80 to scale our black holes 
with charged scalar so that all the solutions have the same value of p. Having fixed all our 
solutions to have the same chemical potential, we now find the temperatures at which two 
of our solutions have the same free energy. A phase boundary appears there if no other 
solutions have a lower free energy at that temperature. We then repeat this with various 
values of /i. The results are shown in Figure El 



^We have dropped the I/IGttG factor in front of the action and set L = 1. Keeping these factors, ah 
of the free energies are proportional to L^ /G which is 0{N'^) in the standard mapping to super Yang-Mills 
theory. This is different from the Hawking-Page transition in global AdS where the free energy is 0(1) at 
low temperatures. 
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Figure 6: Starting from the top left, the phase diagrams for q = 5, q = 2, q = 1.2, q = 1.1, and 
q = 1. In the q = 1 phase diagram, there is a phase boundary between the black hole 
superconductor and the soliton superconductor at about T ~ .001. Here, we have scaled 
the period of the r] coordinate to be 7 = vr. Note that for q = 1.1 and q = 1, there is a 
range of chemical potentials (1.69 < A* < 1.71 for q = 1.1 and 1.74 < /i < 1.81 for q = 1) 
where we find a new superconductor/insulator phase boundary. 

Note that the phase boundary between the non-superconducting phases of the AdS black 
hole and the AdS soliton is known analytically. At zero chemical potential, the boundary 
always occurs at T = l/vr. Also, if there were no superconducting phases, the phase boundary 
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between the AdS black hole and the AdS soliton at zero temperature would be at a chemical 
potential jj, = 2^/^3^''^ ~ 1.86, in our units. 

Unlike the analysis in the probe limit, the superconducting soliton phase seems to exist 
at all chemical potentials. This is different from the case without scalars. A possible physical 
explanation for this may be due to the fact that the curvature of these superconducting black 
hole becomes infinite at zero temperature (analogous to the four dimensional case [20]) and 
the system is trying to avoid singular configurations. We will discuss this further in the next 
section. 

It is clear from Figure |6] that as we lower q, the phase boundaries behave as follows: 
(1) The line dividing the soliton insulator from the soliton superconductor moves to larger 
/i. This is expected since a smaller charge requires a larger chemical potential to cause the 
scalar field to condense. (2) Similarly, the critical temperature for the black hole to develop 
scalar hair is lower as we reduce q. Thus, the slope of the T/fi line marking this phase 
boundary is reduced. (3) The boundary between the black hole superconductor and soliton 
superconductor moves to lower temperature. 

Interestingly, we find a peculiar phase structure for sufficiently low q (e.g. q < 1.1). 
As one lowers g, the two triple points eventually merge into a quadruple pointo and then 
separate again. Once they pass through each other, there is a small range of chemical 
potentials where a black hole conductor can be cooled to form a black hole superconductor 
and then cooled further to become a soliton insulator. That is, there is a new first order 
superconductor/insulator phase transition that opens up at low temperature. The free energy 
as a function of temperature is plotted in Figure [7] at a fixed value of fi in this window. 

6 Summary 

We have studied the four phases of the five dimensional Einstein-Maxwell-scalar theory (12. ip 
corresponding to black holes (with and without scalar hair) and solitons (with and without 
scalar fields). We found the solutions (with backreaction) for various charges q for the scalar 
field. The black holes are always unstable to forming scalar hair along a line of constant T/fi 
exactly analogous to previous studies of the four dimensional case [6] . The soliton is always 
unstable to turning on the scalar field along a line of constant fi as first shown in [11] for 
large q. These phase transitions are usually second order, however we found that the latter 



■^At the quadruple point, the phase boundary between the two soUton phases is second order, which can 
be seen by comparing Figure 2 with Figure 6. 
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Figure 7: The free energy density vs temperature for the q = 1 system at fixed /j, ~ 1.79. The blue 
horizontal line at -F = — 1 is the soliton insulator. The dashed purple line at F ~ —0.97 
is the soliton superconductor. The dotted green line and the thick red line are the black 
hole and black hole superconductor, respectively. Note that as we cool the system from 
high temperature, the free energy favors the black hole, then the superconducting black 
hole, then the soliton insulator. 



becomes first order for small q. 

The complete phase diagrams were constructed for q > 1 and shown in Figure 6. In 
all these cases, the black hole is never the dominant configuration at low temperature. As 
mentioned above, a physical reason for this could be that the curvature of the hairy black 
hole becomes infinite at T = 0, and the system is trying to avoid singular configurations. 
One can test this by including a A|^|^ term in the action (12. ip . With the added term, the 
hairy black hole will be nonsingular at T = 0. For small A, the curvature is still large at low 
temperatures and the phase diagrams are expected be similar to those in Figure 6. But for 
larger values of A, the curvature is smaller at low temperatures and it is not clear whether 
or not the hairy black holes will be preferred at T = 0. 

Although high curvature may be the reason the system with q > 1 prefers the soliton 
over the black hole at low temperature, it seems likely that this will change for g < 1. As 
we discussed in the previous section, the phase boundary, /x = fie, between the soliton and 
soliton superconductor increases as q is lowered. But once /ic > 1.86, the extremal Reissner- 
Nordstrom AdS black hole has lower free energy than the soliton. Since the hairy black hole 
will have even lower free energy, it seems clear that this solution must dominate the phase 
diagram for T = and 1.86 < fi < fic- In this case, the system has to approach the singular 
configuration since no other phase existqj. Higher numerical precision is required to compute 



*An example of a zero temperature phase transition between a soliton and black hole is given in |21) 
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the g < 1 phase diagrams exphcitly. 

Perhaps the most surprising feature of the phase diagrams is the existence of a new 
phase boundary between the hairy black hole and AdS soliton which opens up for q < 1.15. 
This corresponds to a first order phase transition in which a superconductor turns into an 
insulator as the temperature is lowered. It would be interesting to know if any real materials 
behave this way. 

Finally, by considering the case of a neutral scalar field, we constructed the gravitational 
dual of a Bose Einstein condensate of glueballs, and found the maximum energy that can be 
supported this way without thermalizing. 
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